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Abstract 

We discuss the renormalization of local, scalar and pseudoscalar dimension- 
5 operators containing a heavy and a light quark field at scales below the 
heavy-quark mass, using the formalism of the heavy-quark effective theory 
(HQET). We calculate the anomalous dimensions of these operators and 
their mixing to one-loop order. We also perform the one-loop matching of 
gluon and photon penguin operators onto operators of the HQET. We dis¬ 
cuss applications of our results to the mixing of gluon and photon penguin 
operators at low renormalization scales, and to the calculation of l/m-Q 
corrections to meson decays constants. 
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1 Introduction 


Local operators containing both heavy and light quark helds exhibit an interesting 
behaviour under renormalization at scales below the heavy-quark mass mq. Then 
there arise large logarithms of the type ln(mQ/p), characterizing the exchange 
of gluons that are “hard” with respect to the light quark but “soft” with respect 
to the heavy quark. To leading order in an expansion in l/mq, such gluons 
see the heavy quark as a static colour source. The large logarithms can thus 
be summed to all orders in perturbation theory using an effective theory for 
static heavy quarks, the so-called heavy-quark effective theory (HQET) i-i- 
In the HQET, the 4-component heavy-quark held Q{x) is replaced by a velocity- 
dependent 2-component held h^{x) satisfying fh^ = h^, where v is the velocity 
of the hadron containing the heavy quark. Because of the particular hierarchy of 
mass scales involved, the renormalization of heavy-light operators in the HQET is 
called “hybrid” renormalization. Operators in the ehective theory have a diherent 
evolution than in usual QCD. For instance, whereas the vector current q'y^Q is 
conserved in QCD (i.e. its anomalous dimension vanishes), the corresponding 
current q'^^hy in the HQET has a non-trivial anomalous dimension which 

governs the evolution for scales below the heavy-quark mass. 

In the literature, hybrid renormalization has been discussed extensively for 
local current operators of dimension 3 i-i and 4 ||I0| , pd| , as well as for 4- 


quark operators such as the ones governing B-B mixing |^-|^. Here we shall 
consider the renormalization of local dimension-5 operators in the HQET. The 
matrix elements of such operators appear, for instance, at order l/rug in the 
heavy-quark expansion of weak transition form factors. In particular, they con¬ 
tribute to the decay constants of heavy mesons, which have been explored already 
in great detail at leading and next-to-leading order in l/mq |]^. The same op¬ 
erator matrix elements also determine certain moments of meson wave functions 
, which play an important role in the heavy-flavour phenomenology. The the¬ 
oretical predictions for weak decay form factors involve operator matrix elements 
renormalized at the scale mq. Our results can then be used to rewrite these 
matrix elements in terms of ones renormalized at a scale qi <C mq, which may be 
identihed with the scale at which a non-perturbative evaluation of these matrix 
elements is performed (such as the inverse lattice spacing in lattice held theory, 
or the Borel parameter in QCD sum rules). 

Our results also apply to the hybrid renormalization of genuine dimension-5 
operators in QCD, such as the gluon and photon penguin operators, which appear 
in the weak effective Hamiltonian renormalized at the scale mb 0:0 


HeS = Cginib) gss{l + 'y5)o'fiuG^''b + c^(mb) e s (1 -F + ... . (1) 

^Our definition of the coefficients Cg and c.y contains a power of the 6-quark mass, which is 
usually included in the definition of the operators. 
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We will derive the effective Hamiltonian at a lower scale /r < mb, which is rel¬ 
evant to the calculation of hadronic matrix elements of the penguin operators. 
Lattice calculations of such matrix elements, in particular, are usually performed 
in the static theory (HQET), since the b quark is too heavy to be described as 
a dynamical field on present-day lattices. It must be stressed, however, that 
the validity of the effective theory is restricted to the kinematic region where, in 
the rest frame of the heavy hadron, the light quarks and gluons carry momenta 
much smaller than mg. For two-body decays such as B ^ K*'y, one thus needs 
to evaluate the hadronic matrix elements for unphysical particle momenta (i.e., 
|px*| < 0(1 GeV) rather than the physical value |pi<:*| = ~ ^‘k*)/^b — 

2.56 GeV) and then continue the results to the physical region. In more compli¬ 
cated processes such as H —>• > Xg this restriction no longer applies, 

and our results are of direct relevance to the region where the strange particle 
carries a small momentum. 

In Sec. H, we calculate the one-loop anomalous dimension matrix of local 
dimension-5 heavy-light operators carrying zero total momentum. Using the 
equations of motion, the problem is reduced to the mixing of two operators 
containing the gluon field-strength tensor. We solve the renormalization-group 
equation (RGE) for the scale dependence of these operators in the leading loga¬ 
rithmic approximation. In Sec. we extend the basis to the general case where 
the operators carry non-zero momentum. Then there appear four additional op¬ 
erators, which can be written as the total derivatives of some lower-dimensional 
operators. We show that, with a suitable choice of the basis, there is no mix¬ 
ing between these new operators and the ones considered in Sec. |^. In Sec. ||, 
we extend the basis further by including operators containing the photon field, 
and we calculate the mixing between gluon and photon operators under hybrid 
renormalization. This extends the analysis of the mixing of gluon and photon 
penguin operators, which has been discussed previously for scales larger than the 
6-quark mass to low renormalization scales. In Sec. ^ we calculate the 


one-loop matching of the QGD penguin operators onto their HQET counterparts. 
This provides a test of our results for the hybrid anomalous dimensions. In addi¬ 
tion, the exact one-loop expressions for the Wilson coefficient functions may be 
more appropriate to use than the leading-order renormalization-group improved 
results in cases where ln(mQ//i) is not a particularly large parameter; at least, 
they provide an estimate of the importance of next-to-leading corrections. Also, 
the one-loop matching conditions at the scale mg will eventually be part of a full 
next-to-leading order calculation, once the two-loop anomalous dimensions of the 
operators will have been calculated. In Sec. P, we apply our results to the anal¬ 
ysis of higher-order corrections to meson decay constants, and to the discussion 
of moments of meson wave functions. Section contains the conclusions. 
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2 Anomalous dimensions 


We start by considering local, Lorentz-scalar operators of dimension 5, carrying 
zero total momentum. A basis of such operators in the HQET can be con¬ 
structed by considering operators containing two covariant derivatives acting on 
the heavy-quark held. It is convenient to adopt the background-held formalism 
[pl[| , so that it suffices to consider gauge-invariant operators. They are of the 
form where T G We dehne: 


Oi = , 

O2 = gsqifiVyiG^'-'K, 

O 3 = q{iv-Dfhy , 

O4 = qiiDfh^ , 

= q[ipiv-D + iv-Dip]h^ , (2) 

where igsG^’' = is the gluon held-strength tensor, and a^u = 

A basis of pseudoscalar operators 01 can be obtained by replacing q —> ^ 75 . 
As long as we neglect the mass of the light quark and work in a regularization 
scheme with anticommuting 75 , nothing changes by this replacement. Below, we 
shall thus consider scalar operators, but it is understood that the light-quark held 
q could be replaced by 75 ^ or, equivalently, by a left- or right-handed held, qi or 

qn- 

Not all of the operators in (^) are independent when the equations of motion, 
q = 0 and iv-D = 0, are used.|^ They give the relations O 3 = 0, O 5 = —O 2 , 
and 


O 4 = q[[{ip)^''-\9s(T^.uG^^) K 
= —-O 1 + total derivatives. 


{idf{qhP + 2id^{qiD^hP 


(3) 


implying that for zero total momentum all operators can be reduced to the oper¬ 
ators Oi and O 2 . That it is legitimate to use the equations of motion to reduce 
the operator basis has been established in Ref. P^ . 

The scale dependence of the renormalized operators Oi and O 2 is governed by 
a 2 X 2 anomalous dimension matrix, which can be obtained by calculating the 1/e 
poles of the matrix elements of the bare operators in dimensional regularization 
(i.e. in d = 4 — 2e: space-time dimensions). The relevant one-loop diagrams are 
shown in Fig. |^. and Zq are the wave-function renormalization constants for 
the quark helds. In a covariant gauge (a = 1 corresponds to the Feynman gauge, 
a = 0 corresponds to the Landau gauge), they are given by 


rv rv 

Z, = l + C^(3-a)^, Zq = l-Gpa^^. (4) 

^We use the notation iD = id + gsA and (*i2)t = —id + ggA. 
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Figure 1: One-loop diagrams contributing to the calculation of the anomalous 
dimensions of the operators Oi and 02, denoted by the square. Heavy-quark 
propagators in the HQET are drawn as double lines. 


A virtue of the background-held formalism is that the gluon held is not renor¬ 
malized, since ZgZ]^^ = 1 [^. We have calculated the diagrams in Fig. in an 
arbitrary covariant gauge, and with arbitrary momentum assignments (however, 
for zero total momentum). The sum of all diagrams is gauge independent, and 
the dependence on the external momenta combines in such a way that the result 
can be expressed in terms of the matrix elements of the basis operators in (^. 
We hnd 





where N is the number of colours, and Cp = — 1/N) is the eigenvalue of the 

quadratic Casimir operator in the fundamental representation. Eliminating the 
operators O4 and O5 by means of the relations (O4) = and (O5) = —(O 2 ), 

we obtain 



These results dehne a matrix Zij of renormalization constants through = 

ZijiPj). Denoting by z\^'* the coefficient of the 1/e pole in this matrix and using 
the relation m 


lij 


d 

da. 


7(1) 


(7) 
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we obtain the anomalous dimensions appearing in the RGE for the renormalized 
operators: 


Oi R 'JijOj 0 . 


( 8 ) 


At the one-loop order, the result reads 

/ N + Cf 


7 = 




Air 


-2N 


V-fiV + fCF 3N-ICf, 


(9) 


It is remarkable that the eigenvalues of the one-loop anomalous dimension 
matrix are given by irrational numbers (in units of as/Air)\ 


7± = ( 2A^-^ ^ 


/5iV2 


ANCf + 


IQCl 

9 


^-|50± 


'1565 


( 10 ) 


We know of no other case where this happens at the one-loop order. In the 
leading logarithmic approximation, the solution of the RGE (^ is given by 




where 


r+ 0 


0 r_ 


R 


-1 


r± = 


'oisirriQ)' 


7±/2/3o 


(/^) 


( 11 ) 

( 12 ) 


= V 

Here Po = — is the hrst coefficient of the (3 function, and V is the matrix 


3 3 

that diagonalizes the anomalous dimension matrix: 

O' 

^ 47r I 0 7_ , 


(13) 


The explicit form of the evolution matrix is 
/i(r+-fr_)-(iV-|G^)A 

U{mQ,n) = 

\ (-|iV+|G^)A 

where A = (r+ — r_)/( 7 _|_ — 7 _). 


-2NA 

l(r+ + r_) + (A-|G^)A, 


(14) 


3 Complete operator basis 

We now consider the general case where the total momentum carried by the 
operators does not vanish. The purpose of this section is to show that even 
then the operator basis (0i,02) closes under renormalization. The reader not 
concerned about this issue can proceed directly with Sec. |. 
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To obtain a complete basis of local, scalar (or pseudoscalar) dimension-5 oper¬ 
ators, we have to consider, in addition to (^, operators with one or two derivatives 
acting on the light-quark held. These differ from the operators considered so far 
by total derivatives. The number of such operators is reduced signihcantly when 
the equations of motion are used. We hnd that a complete basis contains only 
four new operators in addition to Oi and O 2 , and choose them in the following 
form: 


Ti = id^{qiD^^K) , 

T2 = {iv-df{qK) , 

Ts = {idfiqK), 

T 4 = id^ iv-d (qYK) . (15) 


In momentum space, the total derivative of a local operator corresponds to the 
total external momentum carried by that operator and thus does not affect the 
behaviour under renormalization. Therefore, as far as the calculation of ultravi¬ 
olet divergences is concerned, Ti behaves like a dimension-4 operator, while T 2 , 
Ts and T 4 behave like dimension-3 operators. Under renormalization, operators 
of lower “effective dimension” cannot mix into higher-dimensional operators, but 
the contrary is, in general, not true. As a consequence, the Wilson coefficients of 
the operators Oi and O 2 could, in principle, be modihed by the presence of the 
lower-dimensional operators. 

It follows from this discussion that the anomalous dimension matrix governing 
the mixing of the operators Oi and Tj is of the general form 


f = 



( 16 ) 


where 7 is the 2x2 matrix given (at the one-loop order) in (^. The 4x4 
matrix B describing the mixing of the operators Tj among themselves has been 
calculated in Ref. |^. To all orders in perturbation theory, it has the simple 
form 

/I 
0 
0 


R = 7 hi i + 7“ 


VO 


2 

3 

0 

0 

0 


-1 

0 

0 

0 


0 

0 

0/ 


(17) 


where 7 hi is the anomalous dimension of heavy-light current operators of dimen¬ 
sion 3 [0-0, and 7 “ has been dehned in Ref. O- At the one-loop order, one 
hnds 7 hi = — 7 “ = 7 o(Q:s/ 47 r), where 


10 =-^Cp. (18) 

It remains to hnd the 2x4 matrix A describing the mixing of Oi and O 2 into 
the operators Ti. Since these operators contain total derivatives, the calculation 
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Figure 2: One-loop diagrams contributing to the calculation of the anomalous 
dimension matrix A. 

must be performed with non-zero total momentum. However, a simplihcation is 
that the operators Tj have non-vanishing quark matrix elements at tree level. The 
matrix A can thus be calculated by considering a two-point (rather than three- 
point) vertex function. The relevant one-loop diagrams are shown in Fig. |^. For 
arbitrary quark momenta, we hnd the pole terms: 

^Qbare) = (18/) k) + ..., 

(O/'-®) = Cf^ Uq{p) (V + 2n ■ p/ Uh{v, k) + ... . (19) 

Both matrix elements vanish on-shell, implying that there is no mixing of the 
dimension-5 operators Oi and O 2 into the operators Tj, i.e. 

i = 0. (20) 

In other words, the evolution of the operators Tj is disconnected from that of 
the operators Oj. Moreover, we shall see below that (at least to next-to-leading 
order) the operators Tj are not induced in the matching of QCD operators onto 
HQET operators. This is a virtue of our particular choice of the operators in (|^) 
containing derivatives acting only on the heavy-quark held. Hence, from now on 
the operators Tj can be omitted from our discussion. 


4 Hybrid renormalization of penguins 

We now include in our discussion local operators containing the photon held. This 
is most conveniently done by extending the dehnition of the covariant derivative 
— ipstaA^ — ieA^), so that = iggG^’' + ieF^'', where F^’-' is 

the electromagnetic held. Then the two gluon operators Oi and O 2 in (^ are 
supplemented by their photon counterparts: 

07 = eqa^yF^^^hy , 

O 2 = eq^y/.v^iF^’^hy . ( 21 ) 

The problem is to hnd the mixing of these four operators under renormalization. 
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Figure 3: Additional one-loop diagrams contributing to the calculation of the 
anomalous dimensions in the presence of photon penguin operators, denoted 
by the gray square. 


At leading order in the hne-structure constant a, the gluon operators can mix 
into the photon operators, but not vice versa. Thus, we write the corresponding 
4x4 anomalous dimension matrix in the form 


f = 



( 22 ) 


The 2x2 submatrices X and Y can be obtained from the calculation of the 
diagrams shown in Fig. ^ It is straightforward to see that the values of these 
diagrams can be obtained from the results for the corresponding diagrams in 
Fig. I by performing a simple replacement of colour factors, namely by taking 
the limit N ^ 0 keeping Cp fixed. The first two diagrams in Fig. ^ determine 
the matrix Y and yield 

^ /V rv 

F = 7hil = /(-3Cp)l. (23) 

47r 

Since QCD is blind to the photon held, the operators O 3 and O 4 renormalize 
as dimension-3 heavy-light current operators, i.e. their anomalous dimension is 
given by 7 hi. The submatrix X is obtained from the third diagram in Fig. ^ 
whereas the last diagram vanishes. The result is 

= 1). (24) 

where Qq is the electric charge of the light quark in units of e. Note that the sum 
of these two matrices reproduces the “abelian part” of the matrix 7 in (^), i.e. 

X/Qq + Y = \imN^o 7- 

The three eigenvalues of the anomalous dimension matrix F in units of cts/dvr 
are 7 ± in ([l 0 |) and 70 in ([T 8 |) . The solution of the RGE for the four operators 
Oi can be obtained either by diagonalizing the 4x4 matrix F directly, or by 
solving first the homogeneous equation for the two operators Oj and and 
then inserting the solution for these operators into the inhomogeneous equation 
for the operators Oi and 02- In leading logarithmic order, we hnd that the photon 
operators renormalize multiplicatively: 

OiimQ) =roO]{fx); i = 1,2, (25) 
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with 


7o/2/3o 


To = 


' asirng) ' 


The evolution of the gluon operators is given by 


Oi 

O 2 




rriQ 


Oi 

02 


+ W{mQ,ix)X(0) 


OJ 

01 


(26) 


(27) 


fi \ ^ / fj. 

where X(o) denotes the matrix X in units of a^/dTr, the evolution matrix U has 
been given in (P!2[) and (|^), and 


WimQ,^) = V 


( r+-ro 
7+ -7o 


0 

r_ - ro 


\ 


1 / 


-1 


(28) 


7 


7- - 7o 

Consider now the effective Hamiltonian in (|^. At the scale thq, it contains 
the gluon and photon penguin operators 


Qg = gsq, Q^ = eq , (29) 

with known coefficients Cgirng) and c^^rng). Once again q may be replaced by ql 
or qr. Our goal is to evolve the Hamiltonian down to a low renormalization scale 
fi -C mg. Each of the two penguin operators (renormalized at the scale mg) has 
an expansion in terms of the HQET operators renormalized at a scale fi < mg. 
We dehne a set of coefficient functions by: 

Qg{mg) Oi{mg) = ^ [Cf{fi)Oi{fi)+Cl{fi)Oj] , 

i=l,2 

0'l{mg) = ^ D]{fi)0]{fi). (30) 

i=l,2 

At leading logarithmic order, the initial values of the coefficients at the scale mg 
are determined by a tree-level comparison of operators matrix elements in QCD 
and in the HQET and are thus simply given by Ci{mg) = D'l{mg) = 1; all other 
coefficients vanish at the scale mg. The values of the coefficients at a lower scale 
fi < mg can be read of from the solution of the RGE in (P^D and O- We obtain 

D7(/i) = ro, /^7h)=0, (31) 


and 


cUm) 

c’M 

ciw 


r_|_ -|- r_ 
2 





C^{fi) , 


-2N 


r+ 


r_ 


2QgC 



r_ — rn 

+- 

7- - 7o 

r+ - rp 
7+ - 7o 


+ 2 (^1 

r- - rp 
7- -7o 



C]{fi ), 


(32) 
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With N = 3 and Cp = 4/3, this yields: 

Cf(/i) ~ 0.6966 r_ + 0.3034 r+ , 

Clifx) ~ 0.9652 (r_ - r+), 

C'lifx) ~ Qq (0.7093r_ + 0.0600r+ - 0.7693 tq) , 

Cj(/i) ~ Qq (0.2661 r_ - 0.1355 r+ - 0.1306 tq) . (33) 

As an example, we evaluate the coefficients for the b quark at the scale /i = 1 GeV, 
using Uf = A flavours (i.e. Pq = 25/3), as{mh) = 0.210 and as{^) = 0.458. This 
gives Dj ~ 1 .21, as well as Cf ~ 0.82, Gf ^ 0 .22, CJ ~ —0.26 Qq, and C] ^ 
—0.01 Qq. We observe that the effects of hybrid renormalization are typically of 
order 20%, except for the coefficient C], which is of order ln(mQ//i)]^. 


5 One-loop matching 


In the above discussion of operator evolution, we have combined the tree-level 
matching of the QCD penguin operators onto the HQET penguin operators at the 
scale rrih with the evolution equations solved in the leading logarithmic approxi¬ 
mation. This approach is justihed if ln(mQ//r) S> 1. As an alternative, we shall 
now discuss the full one-loop matching of the QCD operators onto the HQET 
operators. Our calculation will be equivalent to that performed by Eichten and 
Hill for the chromo-magnetic operator containing two heavy-quark helds |^. The 
advantage of this approach is that we will be able to include non-logarithmic 
terms of 0{as). On the other hand, we will not be able to resum logarith¬ 
mic terms to higher orders in perturbation theory. A consistent combination 
of one-loop matching with renormalization-group improvement would require to 
calculate the operator anomalous dimensions to two-loop order, which is beyond 
the scope of our work. However, once this calculation will have been done, the 
one-loop matching conditions at the scale mq derived here will be part of the full 
next-to-leading order analysis. 

To derive expressions for the Wilson coefficients, we must compare the matrix 


elements of the QCD operators Qg and Q.y in (^9]) with the corresponding matrix 
elements of the HQET operators Oi and O]. Order by order in perturbation 
theory, the comparison (“matching”) of these matrix elements dehnes the coeffi¬ 
cient functions. By construction, the differences between matrix elements in the 
two theories are insensitive to any long-distance properties, such as the nature 
of the infrared regulator or of the external states. Therefore, it is legitimate to 
perform the matching calculation using (on-shell) quark and gluon states, and 
working with any infrared regularization scheme that is convenient. Following 
Refs. [|^, 1^, we choose to regulate both ultraviolet and infrared divergences us¬ 
ing dimensional regularization. Moreover, we expand the resulting expressions 
for the Feynman amplitudes in powers of the external momenta, and then set 
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the external momenta to zero inside the loop integrals. Then the only mass 
scale remaining in the QCD calculation is the heavy-quark mass, and hence only 
diagrams containing a heavy-quark propagator in a loop contribute. Likewise, 
in the HQET there is no mass scale left after the external momenta are set to 
zero, and hence all loop integrals vanish on dimensional grounds. So only the 
tree-level matrix elements of the HQET operators multiplied by their Wilson co¬ 
efficient functions remain. That is why this particular regularization scheme is 
most economic for our purpose. 

In the matching calculation, we have to consider all operators in the HQET 
that have the same quantum numbers as the QCD operator in (|^) . In addition 
to the genuine dimension-5 operators Oi and O] and the total derivatives Ti 
encountered so far, this set includes in principle also operators proportional to 
the heavy-quark mass mq. They are: 


Si = mqiv-d {qK ), 

5*2 = mQid^{qYhv) , 

S'3 = mqqiv-Dhv, 

S '4 = mqqiph ^, 

S 5 = mqqhy. (34) 

The matrix elements of the operators S 3 and S '4 — S 2 evaluated between physical 
states vanish by the equations of motion. Nevertheless, these operators may be 
induced when the matching calculation is performed with unphysical states, such 
as on-shell quarks and gluons. 


/i, a 



mv k p 


Figure 4: One-loop diagrams contributing to the matching calculation for the 
gluon and photon penguin operators, Qg and Q^, denoted by the white and 
gray squares. 


Since some operators (T2, T3 and T4, as well as Si, S 2 and S' 5 ) do not have 
gluon matrix elements at tree level, we have to perform the matching calculation 
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by considering both the qnark-qnark and the qnark-qnark-glnon vertex fnnc- 
tions, with arbitrary (bnt on-shell) external momenta. The relevant diagrams 
are shown in Fig. We have performed the calcnlation of these diagrams in an 
arbitrary covariant gange, and hnd that the resnlts for the coefficient fnnctions 
are gange independent. The on-shell wave-fnnction renormalization constant for 
the heavy-qnark held in onr regnlarization scheme is 




ml qrf-l)r(2-rf/2) 

Anfj?) {d — 3) 


(35) 


which is also gange independent. The light-qnark held is not renormalized in this 
scheme (i.e. Zq = 1). Note that only the second diagram in Fig. ^involves a light- 
qnark propagator in a loop. It is easily seen that the contribntion of this diagram 
is linear in the external glnon momentnm, bnt independent of the momentnm 
of the light qnark. Hence, it follows that no contribntions proportional to the 
light-qnark momentnm appear in the matching calculation. This means that the 
operators T, and Si, S 2 are not induced by radiative corrections, at least not to 
0{as). 

Let us now present the result of the evaluation of the diagrams involving the 
gluon penguin operator Qg. For the sum of all contributions, we hnd 


(Q,) = Cf (O,) + Cl {O 2 ) + Cl {01) + Cl {CD 

+ Cl (Os) + Cs (-2{Ss) + (S4> + 2 (s,)) , (36) 


where the exact expressions for the coefficients in an arbitrary space-time dimen¬ 
sion d are 


Cf 

d 

d 

d 


1 + y [{2d^ - 23d^ + 71d - 62 )Cf - 


-K 


2{d - 2){d - 3){d - 4 :)Cf - d 


KQqCF {d — 3){dS — Qd + 16), 
-2KQqCF {d - 4)(d2 -Qd+ 10), 


d - lOd^ + 25d - U)N 
- 9d‘^ + 23d - 14)Arj , 


(37) 


as well as 


d = -4 :KCf {d - l)(d - 4), 

Cs = -2 KCf {d - l){d - 3), 

where we have abbreviated 

_ ^ (luLd'' r(2 - d/2) 

dvr ydTT/i^y {d — 2){d — 3) 


(38) 


(39) 
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Note that the operators O 3 and S 3 induced by one-loop matching in (^) 
vanish by the equations of motion. These operators appear only because the 
matching calculation is (legitimately) being performed using unphysical quark 
and gluon states, rather than physical hadron states. The quark states satisfy the 
on-shell conditions v ■ k = 0 and ^ = p'^ = 0 , which differ from the true equations 
of motion at 0{gs)- The presence of operators that vanish by the equations 
of motion in (|3^ is, however, without physical signihcance. For all practical 
applications these operators can be ignored, since their physical matrix elements 
vanish. A similar statement obtains for the operators 5*4 and S 5 . That fact that 


they appear in (|36|) in the combination S '5 -|- 154 is a consequence of the so- 


called reparametrization invariance of the HQET [^. Indeed, this combination 
is nothing but the HQET counterpart of the QCD operator m^qQ. But this 
operator has the form of an off-diagonal mass term, which has no observable effect 
as it can be removed by a redehnition of the quark helds. This simply shifts the 
physical quark masses by an amount of OiCpasm^). This discussion shows a 
posteriori that, as shown in (|5(]|), the matching of the QCD operator Qg onto the 
HQET operators involves only genuine dimension-5 operators, preserving thus 
the structure of an effective held theory, in which operators enhanced by powers 
of the large scale mq do not contribute to physical matrix elements. 

A similar calculation for the matrix element of the photon penguin operator 
gives, omitting now unphysical operators. 


(S,) = B? {OJ} + BJ (OJ), 


(40) 


where 


DJ = 1 + j Cf {d^ - I5d + 34), 
= -2KCf {d - 4)2 . 


(41) 


Let us now set d = 4 — 2e, take the limit e —0 and remove the poles 
in 1/e using the Ms subtraction prescription. In doing this, we have to take 
account of the fact that not only the matrix elements of the HQET operators 
are ultraviolet divergent, but also the matrix elements of the QCD operators Qg 
and themselves. If we are interested in the evolution of the operators below 
the scale mq, taking their values at /i = mq as given, we have to remove this 
second type of divergence by renormalizing the bare QCD operators at the scale 
/i = mq. This is accomplished by a 2 x 2 matrix which is given by |^, p!9|:P| 


z-^ = 1 + 


cr c 


m: 


Q 


Stt \ dvr/i 


d/2-2 / 2 Cf 

2 , r (2 - d/ 2 ) 


0 


(42) 


\SQqCF IOCf — 4A^ ^ 


^The matrix in parenthesis contains the one-loop anomalous dimensions of the QCD penguin 
operators, after a factor of the heavy-quark mass has been removed. 
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which multiplies the 2x4 matrix of coefficient functions 


/ /^9 r-^9 
Oi L/2 O 2 


c = 


( 43 ) 


VO 0 DJ D] 


from the left. The remaining divergences are removed by subtracting the Xje 
poles using the Ms scheme. The result is: 


and 


D7(/i) = \-Cp 


a.. 



^I(h)=0, 


(44) 


Cf(/i) = 1 + (AT + C^) ^ In^ + (iV - ^ , 

OTT TUq \ 4 / TT 

C|(h) = ln-4-(iV-2C'i.)^, 


47r 


m 


Q 


271 




CIM = Q,Cf 


(y.a 


TT 


(45) 


The Wilson coefficients obey the RGE 


d/i 


(46) 


where T has been given in (^^-(^). The fact that our explicit results for the 
coefficient functions satisfy this equation is a strong check of our results. 

As at the end of Sec. we now evaluate as an example the coefficients for the b 
quark at the scale = I GeV, using rrib = 4.8 GeV and Ug = ag^y/JTmh) = 0.282. 
This gives Dj ~ 1.02, as well as Cf ~ 0.97, Gf ^ 0.20, G7 ^ —0.22Qq, and 
C] — 0.12 Qq. The comparison with the leading-order renormalization-group 


improved results presented after (^^ gives an idea of the importance of next-to- 
leading corrections. In some cases, such as Dj and C 2 , the non-logarithmic terms 


of 0{as) in (^) are quite important. For those coefficients a full next-to-leading 
order calculation would be desirable. 


6 Meson decay constants and wave functions 

Another application of our results is the calculation of higher-order corrections in 
the heavy-quark expansion of current matrix elements. As an example, we discuss 
the matrix elements of heavy-light currents between a ground-state meson and 
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the vacuum, which dehne the meson decay constants Jm- The heavy-quark ex¬ 
pansion for heavy-meson decay constants has been discussed at next-to-leading 
order in l/mq in Ref. [^. In general, at order l/rriq in the expansion there 
appear contributions from the matrix elements of local current operators of di¬ 
mension 3 -|- n, and (for n > 1) from the matrix elements of non-local operators 
containing the time-ordered products of lower-dimensional current operators with 
terms from the effective Lagrangian of the HQET. Our results are relevant to the 
discussion of the local corrections of order 1/rriq. They can be expressed in terms 
of the matrix elements of operators of the typeQ qT iD^iD^hy. Unlike the case 
considered so far in this paper, these operators transform as vectors and axial 
vectors under the Lorentz group. However, our results still apply because of 
heavy-quark spin symmetry. In fact, the most general matrix element of opera¬ 
tors of the type shown above can be written in the form 


(0|grzDWh,|M(u))^ = -vWF(/i)Tr{0^"r>l(u)} , (47) 


where F ^ /m^/Fim is the leading contribution to the meson decay constant, 
and 

.... 1 ^ f —z 75 : pseudoscalar meson M(v), . 

= ^ ^ (48) 

2 ; vector meson M (e, u), 

is a Dirac matrix representing the spin wave function of the ground-state mesons 
in the HQET [^. Here v is the meson velocity, and e is the polarization vector 
of the vector meson {e ■ v = 0). The matrix M.[v) satishes = —A4{v). 

The most general decomposition of the tensor form factor 0^^ consistent with 
Lorentz covariance, heavy-quark symmetry and the equations of motion involves 
the “binding energy” A = rriM—fnq, which is a scale-independent mass parameter 


of the HQET, as well as two new parameters and Xjj |TB|: 


0 /^^ = 2 


A^ + Ai(/i)+A|(/z) {v>^v'^-gn 


+ 2A|(/i) + vn + ^ A^(p) (1 - ^). 


(49) 


Heavy-quark spin symmetry ensures that the relation (^^ holds for any matrix 
T. Thus, we are free to set T = 75 ( 7 ^^^ — 'yu'Vfj.) or T = in which case we 

End relations between the parameters A|; and A^ and the matrix elements of the 
pseudoscalar operators and evaluated between a pseudoscalar meson 
and the vacuum: 


( 0 I \M{v)) = -i^/rn^ F(p)A|(p), 

( 0 I ioP - of \M{v)) = -zV^F(p)Ai(p). (50) 

^There are also contributions from operators that are total derivatives; however, their matrix 
elements can be expressed in terms of parameters already encountered at order 1 /toq in the 
expansion. 
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In the rest frame of the hadron, these are matrix elements of pnrely chromo¬ 
electric and chromo-magnetic operators, respectively. 

The scale dependence of these matrix elements is determined by the scale 
dependence of the operators Oi and O 2 , which is described by the evolution 
matrix given in We obtain 


ro AKmg) = 


r I -|- r_ 


JY r 


r_ 




7+ -7- 


(/^) 


N 


C, 


r: 


r_ 




7+ 

2 


7- 


H 


(/^) 


~ (0.7413r_ + 0.2587r+) + 0.4379 (r_ - r+) A^(/x), 


ro^nirng) = 



W 

2 


-Cf 

3 


r+ — r_ 


A^ 


7+ - 7 - 

~ (0.7413 r+ + 0.2587r_) A^(/x) + 0.4379 (r_ - r+) A|(/x), 


e(/^) 

(51) 


where the factor tq on the left-hand side compensates for the scale dependence of 
the parameter F(/i), according to the equation F^tiiq) = roF(/i), with vq given 

Let us now discuss the specihc example of meson decay constants to leading 
logarithmic order. The local corrections in the heavy-quark expansion are ob¬ 
tained from the expansion of the flavour-changing current ^7“(1 — 75 )Q, using 

i 

ip — iv-D j gs 


Q —-\- 


2m 


Q 


+ 5 —2 iG^^K + ... . 


4>"q 


(52) 


The corresponding expressions for the meson decay constants are 


\ l-d 


A 


M 


6m, 


+ 


Xl 


mg) 


Q 


12 ml 


+ 


-|- non-local terms 


(53) 

where the spin-dependent coefficient cIm takes the values 3 and —1 for pseu¬ 
doscalar and vector mesons, respectively. The non-local matrix elements of order 
1/mg are discussed in Ref. [^. The new ingredient in (|5BD is the local correction 
of order l/mg, which is determined by the single parameter X\ renormalized at 
the scale mg. Using our result (^l|), this parameter can be related to the val¬ 


ues of X\ and A|^ renormalized at some lower scale /x, which is typically the 
scale intrinsic to some non-perturbative calculation of these hadronic parame¬ 
ters. For instance, QCD sum rules have been used to predict that A|;(/i) ~ 
0.11 GeV^ and A|^(/x) 0.18 GeV^ at the scale /x 1 GeV [|T^. For the b 
quark, e.g., using the numbers given at the end of Sec. this translates to 
X\{mb) ~ 0.69A|;(/x) -|- 0.08A^(/i) ~ 0.09 GeV^, which is the value to be used 
in (^). The corresponding correction to the R-meson decay constants is clearly 
very small, indicating a good convergence of the local terms in the heavy-quark 
expansion. Before a meaningful prediction for the decay constants at order l/m/g 
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can be obtained, it wonld however be necessary to inclnde the non-local terms as 

well. 

For completeness, we note that the hadronic parameters and \\ are of a 
more general interest, since they are related to the second moments of meson wave 
fnnctions (y9±(a;,p). Introdncing a vector 2 ; on the light-cone v ■ z = t) 

and working in light-cone gange (^4+ = 0), we dehne 


^ J dfe*""* (0 I g(0)7±r h^{z) |M(n))^ = ij,)Tt {'y±T M{v)} , 

(54) 

where 7 ± are the light-cone projections of the Dirac matrices. Dehning the mo¬ 
ments 


(w">± = jdul ip±(ul, fi) ul ", 


(56) 


which are normalized snch that {uj^)± = 1, one can show that ||l6| p| 


{‘^)± = ^ 5) 

“ (5 * 5) + (5 + 5 


Onr resnlts in (|^ can be used to control the scale 
moments. 


dependence of the second 


7 Conclusions 

We have discussed the hybrid renormalization of local, scalar and pseudoscalar 
dimension-5 operators containing a heavy and a light quark held. Our results 
determine the scale dependence and mixing of such operators at scales below the 
heavy-quark mass. We have calculated the corresponding anomalous dimensions 
at the one-loop order, which is sufficient to obtain the renormalization-group 
evolution of the operators in the leading logarithmic approximation. 

Two applications of our results have been discussed in detail. The hrst one 
concerns the renormalization of genuine dimension-5 QCD operators at low renor¬ 
malization scales. Important examples are the gluon and photon penguin oper¬ 
ators, which appear in the effective weak Hamiltonian renormalized at the scale 
rub, and whose evolution at high scales /r S> mb is well known. We have dis¬ 
cussed the mixing of these penguin operators at a low scale n <^mb and derived 
the effective Hamiltonian for this case. This is relevant for the calculation of 
hadronic matrix elements of the penguin operators performed, e.g., using lattice 

^These results are valid in a regularization scheme without a dimensionful regulator, such 
as dimensional regularization. 
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simulations. In addition to the evolution at leading logarithmic order, we have 
also performed the full one-loop matching of gluon and photon penguin operators 
onto their effective-theory counterparts. Besides providing a test of our results 
for the operator anomalous dimensions, this calculation will eventually be part 
of a full next-to-leading order renormalization-group improved analysis, once the 
two-loop anomalous dimensions of the operators will have been calculated. 

The second application concerns the calculation of certain higher-order correc¬ 
tions in the heavy-quark expansion of current matrix elements, such as they ap¬ 
pear in the description of weak decay form factors. In particular, local dimension- 
5 operators appear at order l/uig in the heavy-quark expansion of meson decay 
constants. We have discussed, as an example, the local l/rug corrections to the 
decay constants of heavy mesons at the scale mg, at which a single new hadronic 
parameter appears. Our results can then be used to evolve the result down to 
a lower scale p, at which non-perturbative evaluations of the relevant hadronic 
matrix elements may be performed. A similar analysis could be performed for 
semileptonic transition form factors describing, e.g., semileptonic decays such as 
B Triu. 
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